Topological phase transitions for interacting finite systems 
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In this paper, we investigate signatures of topological phase transitions in interacting systems. 
We show that the key signature is the existence of a topologically protected level crossing, which is 
robust and sharply defines the topological transition, even in finite-size systems. Spatial symmetries 
are argued to play a fundamental role in the selection of the boundary conditions to be used to 
locate topological transitions in finite systems. We discuss the theoretical implications of this result, 
and utilize exact diagonalization to demonstrate its manifestations in the Haldane-Fermi-Hubbard 
model. Our findings provide an efficient way to detect topological transitions in experiments and in 
numerical calculations that cannot access the ground-state wave function. 

PACS numbers: 03.65.Vf, 71.10.Fd, 05.30.Fk, 05.30.Rt, 



The discovery of the integer quantum Hall effect 1 pro- 
vided a new type of quantum phase transition, the topo- 
logical transition, which does not depend on spontaneous 
symmetry breaking. In recent years, the interest in topo- 
logical states of matter and topological transitions was 
renewed by the discovery of a class of topological insula- 
tors with time- reversal symmetry 2 ' 3 . For noninteracting 
insulators, it is believed that all possible topologically 
ordered states have been classified 4 ' 5 , and the transition 
between topologically distinct states must feature a clos- 
ing of the single-particle gap 2,3,6 . At present, much of 
the recent development has focused on the role of inter- 
actions in topological insulators, specifically the search 
for interaction- induced topological insulators 7-11 and un- 
derstanding the nature of quantum phase transitions be- 
tween different topological classes 12-21 . 

Some of the most important challenges in the study of 
interacting topological insulators lie in developing a clas- 
sification scheme and in the difficulty to accurately com- 
pute the corresponding topological indices in an inter- 
acting system. For a Chern insulator, the standard way 
of calculating the topological invariant (the Chern num- 
ber) involves determining the ground-state wave func- 
tion with twisted boundary conditions 22,23 or by taking 
a three-dimensional (3D) integral of the Green's func- 
tion 24 ~ 26 . Alternatively, it has been proposed that topo- 
logical order can be ascertained by examining the single- 
particle gap near the edge states 13 or the entanglement 



spectrum 



27,28 



Unfortunately, these approaches are com- 



putationally very challenging. 

In this Rapid Communication, we show that, for 
all interacting topological insulators that can be classi- 
fied by the Chern number, the topological transition is 
sharply defined in finite-size systems. In contrast to non- 
interacting infinite systems, where the topological transi- 
tion is always marked by the closing of the single-particle 
gap, we show that, for an interacting finite-size system, 
the interaction-driven topological transition may be char- 
acterized by the closing of the excitation gap without 



necessarily closing the single-particle gap. This closing 
of the excitation gap can be viewed as a topologically 
protected level crossing. Next, we show that for models 
with inversion symmetry the level crossing can only occur 
for boundary conditions that are invariant under inver- 
sion. Based on this finding, we provide a prescription for 
efficiently determining the topological transition and, us- 
ing exact diagonalization, perform a representative study 
of interacting topological states in the simplest case of 
a lattice quantum Hall state with broken time-reversal 
symmetry. 

We begin by formulating a general theory that cap- 
tures such transitions. In contrast to ordinary quantum 
phase transitions 29 , a topological phase transition does 
not require spontaneous symmetry breaking and can be 
precisely defined and observed even in finite-size systems. 
Consider a two-dimensional insulator with Hamiltonian 
H(X), where A is a control parameter. Given twisted- 
boundary conditions 30,31 , the Chern number can be de- 
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{(0^*18^)- (d* y **\d^)), (1) 



where is the exact many-particle wave function and 
4>x {4>y) are twists along the x (y) direction. As long as 
a unique ground state is found for all twisted boundary 
conditions, the integral of Eq. (1) is quantized to an inte- 
ger value for any size system. In other words, regardless 
of the size of the system, we can always define a topolog- 
ical transition between insulators as the place where C 
changes its value from one integer to another. 

Now we adiabatically vary A from Ai to A2. If the 
excitation gap remains finite during this procedure for 
all twisted-boundary conditions, the value of the Chern 
number must remain invariant because the topological 
index is quantized to integer numbers for gapped sys- 
tems. This observation immediately implies that if the 
topological index changes its value, then the excita- 
tion gap Acx = Ei — Eq, with Eq (Ei) the energy of 
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the ground (first-excited) state, must vanish for some 
twisted-boundary condition at the topological transition. 
In direct contrast to an ordinary quantum phase transi- 
tion, where finite-size effects in general result in a finite- 
size gap, this phenomenon of a vanishing excitation gap 
remains even in finite-size systems, where the vanishing 
of Aex implies the existence of a level crossing between 
the lowest two states. We emphasize here that this level 
crossing is required by the topological properties of the 
ground-state wave functions, and thus we refer it to as a 
topologically protected level crossing. 

Identifying this topologically protected level crossing 
point, in principle, requires the computation of excita- 
tion gaps for every twisted-boundary condition. This 
difficulty can be avoided if we focus on a special class of 
topological transitions where (a) the system has space- 
inversion symmetry and (b) the topological index changes 
by an odd number at the transition. These two conditions 
are satisfied in a large class of topological transitions 
[including the Haldane-Fermi-Hubbard (HFH) Hamilto- 
nian 13 , which we investigate below as a test model]. With 
space-inversion symmetry, the excitation gap at ((j> x , 4> y ) 
must coincide with its partner (—<j) x ,—(j)y). At the topo- 
logical transition point, this symmetry relation implies 
that if the gap closes at some (4> Xl <fi v ), then (—<fi Xl —4>y) 
also has a level crossing. Consequently, there are in gen- 
eral an even number of level crossings and the Chern 
number must change by an even integer. In order for the 
Chern number to change by an odd value at the transition 
point, the level crossing must occur at one of the bound- 
ary conditions which are their own space-inversion part- 
ners: (0,0), (tt, 0), (0,7r), and (tt, it). Thus, we only need 
to examine the excitation gap for these high-symmetry 
boundary conditions to identify the topological transi- 
tion. In addition, systems with higher rotational sym- 
metry can simplify this further 31 . 

Finally, we emphasize that the existence of a level 
crossing is a necessary condition for a topological tran- 
sition, instead of a sufficient one. By naively looking at 
the excitation gap, one cannot distinguish the topolog- 
ical phase transition from an accidental level crossing. 
However, if one knows that the topological index does 
change, e.g., by calculating the Chern number or by the 
use of limiting arguments, then the level crossing must 
be associated with the topological transition. 

To demonstrate the physics described above, we use 
a thick-restart Lanczos algorithm 32 to study a model 
which has a lattice quantum Hall state with broken time- 
reversal symmetry 6 . Here, we consider the HFH Hamil- 
tonian 13 , 

H = -t 1 5^( C t Cj .+H.c) 

(2) 

-t 2 J2 (e^cjej + H.c.) +V y £n i n ji 

iij)) {ij) 

on a honeycomb lattice at half-filling, where cj (c { ) are 
the fermion creation (annihilation) operators at site i and 



m = cjcj is the corresponding number operator. Here 
t\ (£2) are the nearest-neighbor (next-nearest-neighbor) 
hopping amplitudes, V is a repulsive nearest-neighbor 
interaction, and the next-nearest-neighbor hopping term 
has a complex phase 4>ij = i</> corresponding to loops 
in the anti-clockwise (clockwise) direction. In what fol- 
lows, we restrict our study to clusters whose symmetry 
in momentum space contains the zone corner k = A', as 
justified in Ref. 13. Also note that we set the unit of 
energy t\ = 1 and the definitions of all the observables 
are presented in the Supplemental Material 31 . 

For small V the system is a gapped topological insu- 
lator, which has a unique ground state and a finite gap 
Aex . In the limit V — > 00, the system turns into a topo- 
logically trivial charge-density-wave (CDW) insulator in 
which all of the particles are located on one sublattice. 
Here, the system has a doubly degenerate ground state, 

i.e., Aox^ = 0, and two finite gaps: the excitation gap 
(2) 

A ex = E2 — Eq, where E2 is the energy of the second 
excited state, and single-particle gap A sp , which is the 
energy required to add or remove a particle from the 
system 31 . In general, the onset of CDW order and the 
change in the topological index may occur at different 
interaction strengths, Vc and Vr, respectively, opening 
up a topological Mott-insulating region. For this Hamil- 
tonian, we observe two cases that are related to the sym- 
metry of the cluster: (I) Vc < Vr and (2) Vc = Vr- 

Figure I depicts the properties of the HFH Hamilto- 
nian for the 24C cluster [see the inset in Fig. 1(d)] with 
parameters that typify the case Vc <Vr- In Fig. 1(a), we 
show the CDW structure factor 31 and the Chern num- 
ber. Here the jump in the structure factor marks the 
CDW transition at V c = 4.022 ± 0.001. In addition to 
the CDW transition, the topological index also changes 
its value as V increases, and we identify the topological 
transition at V T = 4.5281 ± 0.0001. 

Next, we show the four lowest-energy states in 
Fig. 1(b), with an inset focusing on the avoided level 
crossing at V = Vc- In addition, there is a topologically 
protected level crossing at V = Vr (not visible). This 
can be seen more clearly in Fig. 1(c), where we show 

the single-particle gap A sp and excitation gaps A ex and 

(2) (2) 
Aex . Here A sp and A ex both have a pronounced mini- 
mum at V — Vc, where both gaps are expected to vanish 
in the thermodynamic limit. The topological transition, 
on the other hand, is characterized by a vanishing exci- 
tation gap Aex , not necessarily by the vanishing of the 
single-particle gap. This is in direct contrast to Ref. 18, 
which claims that the topological transition is connected 
to the minimum of the single-particle gap. In addition, 
we emphasize that for all of the clusters we studied, the 
closing of the excitation gap at the topological transition 
always takes place for the periodic boundary condition 
case 4>x = 4> v = (Ref. 31). 

One natural consequence of the level crossing seen in 
Fig. 1(c) can be observed in the fidelity metric 31 g(V, SV), 
which has been shown to be a sensitive indicator of quan- 
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FIG. 1: (Color online) (a) CDW structure factor Scdw an d 
Chern number C, (b) ground-state energy Eq and first three 
excited state energies (Ei, E2, and E3) with an inset that 
shows a close up view of the avoided level crossing at V = Va , 

(c) excitation gaps Aex and a single-particle gap A sp , and 

(d) fidelity metric g(V, SV) with SV = 10~ 4 as a function 
of the interaction strength V/t\ for the 24C cluster [see the 
inset of (d)], and parameters t\ = 1.0, t% = 0.8, and <j) = 7r/2. 
Note that this cluster does not possess all of the symmetries 
present in the infinite system. 



turn phase transitions 13 ' 33-36 . We illustrate this quantity 
in Fig. 1(d). While the CDW transition is marked by a 
peak with finite width (independent of SV) indicative of a 
traditional (first-order) phase transition, the topological 
transition is characterized by a singular point where the 
overlap goes to zero, and the fidelity metric has a singular 
peak with height 2/N(SV) 2 and width ~ SV, where N 
is the number of sites. Here we emphasize that because 
there is a topologically protected level crossing, this sin- 
gular behavior of the fidelity metric will always occur for 
a topological transition. Due to the singular nature of 
this peak, numerical calculations of g may easily miss 
this feature, instead observing a jump discontinuity 13 . 

The second case, Vc = Vr, is more representative of 
the model in the infinite limit. Indeed, for any cluster 



FIG. 2: (Color online) (a) CDW structure factor Scdw and 
Chern number C, (b) ground-state energy Eo and first three 
excited state energies (Ei, E2, and E3) with an inset that 
shows a close up view of the level crossing at V = Vc = Vr, 

(c) excitation gaps A ex and a single-particle gap A sp , and 

(d) fidelity metric as a function of the interaction strength 
V/ti for the 24D cluster [see the inset of (d)] and the same 
parameters as Fig. 1. 



which preserves the full symmetry of the honeycomb lat- 
tice we find that Vc = Vr for all parameters studied. In 
Fig. 2, we show the properties of the 24D cluster [see the 
inset in Fig. 2(d)] with the same parameters as in Fig. 1. 
In Fig. 2(a), the CDW transition, marked by the jump in 
the structure factor, and the change in the Chern number 
both occur at V = V C = V T = 3.9813 ± 0.0001. As a re- 
sult, there is a level crossing [shown in Fig. 2(b)] which is 
topologically protected and the ground state is triply de- 
generate. This three-fold degeneracy becomes evident by 

examining the excitation gaps Aex and A^ [Fig. 2(c)], 
which both approach zero as V — > Vr- In the fidelity 
metric [Fig. 2(d)], we see only the singular peak associ- 
ated with the topological transition. In general, we find 
that the two features do coexist, with the CDW fidelity 
peak becoming much sharper and, in many cases, com- 
pletely obscured by the singular peak at the topological 
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FIG. 3: (Color online) tz-V phase diagram for the (a) 12^4, 
(b) 18C, (c) 24C, and (d) 24D clusters with parameters 
ti — 1.0 and 4> = 7r/2. The red squares indicate the onset 
of CDW order, and the blue circles indicate the topological 
transition. The shaded region marks the coexistence of CDW 
and topological order. All points to the vertical line are domi- 
nated by strong finite-size effects. Illustrations of the 12A and 
18C clusters are shown as insets in the corresponding panel. 



transition. 

In Fig. 3, we show the t 2 -V phase diagrams for clusters 
from 12 to 24 sites. In general, the system is a topolog- 
ical insulator (TI) at weak-coupling and a topologically 
trivial CDW insulator at strong-coupling. For clusters 
without six-fold rotational symmetry [Figs. 3(a)-3(c)], we 
find a region of parameter space with a coexistence of 
topological and CDW order. For clusters with six-fold 
rotational symmetry [Fig. 3(d)], we find that the CDW 
and topological transitions always coincide. As this sym- 
metry is present in the thermodynamic limit and given 
the reduction we see in the region where CDW and topo- 
logical order coexist as we increase the size of the clusters 
without that symmetry, we conclude that the topological 
CDW insulator phase does not exist in this model in the 
thermodynamic limit. In addition, we find no evidence 
for a topologically trivial insulating phase without CDW 



order 18 . 

It is known that, at t 2 = 0, the system exhibits a semi- 
metal-CDW transition at finite V in the thermodynamic 
limit. However, in Fig. 3, we observe a sudden decrease 
in Vc for small t$, resulting in a very small value of Vc 
at ti = 0. This apparent contradiction is due to finite- 
size effects, which become dominant as V, t% —> 0. In 
this region, the single-particle gap becomes very small 
(~ t%), so much larger systems (with sizes larger than the 
inverse gap) need to be studied to accurately determine 
the phase boundaries. To highlight this observation, we 
mark the onset of strong finite-size effects (t 2 = L- 1 ) by 
a dashed line in Fig. 3, where L is the linear size of the 
system. The drop of Vc takes place in the regime with 
strong finite-size effects (t 2 < L^ 1 ). As the system size 
increases, the region with t 2 < L~ x is pushed to t 2 = 0, 
so the sudden drop of Vc at small t 2 disappears in the 
thermodynamic limit, resulting in a finite Vc in the limit 

*2 -> 0. 

In summary, we have shown that the closing of the 
excitation gap is a signature of the topological tran- 
sition in interacting systems. This topologically pro- 
tected level crossing exists even in finite-size systems, 
and the resulting singular behavior at the transition can 
be observed in quantities such as the fidelity metric. 
When coupled with the use of spatial symmetries to sim- 
plify the choice of boundary conditions, this phenomenon 
provides an efficient scheme for locating the topolog- 
ical transition whichcan be straightforwardly general- 
ized to time-reversal invariant topological insulators and 
fractional topological states. Aside from a few special 
cases 37,38 , our findings provide a generic methodology 
to locate a topological transition in interacting finite-size 
systems via an experimentally measurable quantity. Con- 
sequently, this allows for the study of topological phases 
and transitions in cold-atom experiments and computa- 
tional approaches in which the ground-state wave func- 
tion is not accessible. 
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